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Abstract 



An analysis of nuclear properties based on a relativistic energy functional 
containing Dirac nucleons and classical scalar and vector meson fields is dis- 
cussed. Density functional theory implies that this energy functional can 
include many-body effects that go beyond the simple Hartree approximation. 
Using basic ideas from effective field theory, a systematic truncation scheme 
is developed for the energy functional, which is based on an expansion in 
powers of the meson fields and their gradients. The utility of this approach 
relies on the observation that the large scalar and vector fields in nuclei are 
small enough compared to the nucleon mass to provide useful expansion pa- 
rameters, yet large enough that exchange and correlation corrections to the 
fields can be treated as minor perturbations. Field equations for nuclei and 
nuclear matter are obtained by extremizing the energy functional with re- 
spect to the field variables, and inversion of these field equations allows one 
to express the unknown coefficients in the energy functional directly in terms 
of nuclear matter properties near equilibrium. This allows for a systematic 
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and complete study of the parameter space, so that parameter sets that ac- 
curately reproduce nuclear observables can be found, and models that fail to 
reproduce nuclear properties can be excluded. 

Chiral models are analyzed by considering specific lagrangians that real- 
ize the spontaneously broken chiral symmetry of QCD in different ways and 
by studying them at the Hartree level. The resulting energy functionals are 
special cases of the general functional considered earlier. Models that include 
a light scalar meson playing a dual role as the chiral partner of the pion 
and the mediator of the intermediate-range nucleon-nucleon interaction, and 
which include a "Mexican-hat" potential, fail to reproduce basic ground-state 
properties of nuclei at the Hartree level. In contrast, chiral models with a 
nonlinear realization of the symmetry are shown to contain the full flexibil- 
ity inherent in the general energy functional and can therefore successfully 
describe nuclei. 
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I. INTRODUCTION 



There is a long history of attempts to unite relativistic mean-field phenomenology with 
manifest chiral symmetry In particular, it has been tempting to build upon the linear sigma 
model fl|-|3[ , where a light scalar meson plays a dual role as the chiral partner of the pion and 
the mediator of the intermediate-range nucleon-nucleon (NN) attraction [f§-[16|]. In earlier 



work [17]], we surveyed a broad class of these chiral hadronic models and observed that they 
fail to reproduce basic properties of finite nuclei at the Hartree level. In this paper, we 
demonstrate the generic failure of this type of model using a more complete approach |18 



for analyzing the relationship between model parameters and nuclear observables. We also 
illustrate the characteristics of chiral models that can successfully describe finite nuclei. 

We base our analysis on an energy functional, which depends on valence-nucleon Dirac 
wave functions and classical scalar and vector meson fields. Extremizing the functional leads 
to coupled equations for finite nuclei and nuclear matter [see Eqs. (fj)-®]. The successes of 
relativistic mean-field models have shown that these variables (nucleon, scalar, vector) allow 
an efficient and natural description of bulk and single-particle nuclear properties |T9|-p7| . 

Although the energy functional contains classical meson fields, this framework can ac- 
commodate physics beyond the simple Hartree (or one-baryon-loop) approximation. This is 
achieved by combining aspects of both density functional theory (DFT) [^-|3T| and effec- 
tive field theory (EFT) |32| p4|. In a DFT formulation of the relativistic nuclear many-body 
problem, the central object is an energy functional of scalar and vector densities (or more 
generally, vector four-currents). Extremization of the functional gives rise to Dirac equa- 
tions for occupied orbitals with local scalar and vector potentials, not only in the Hartree 
approximation, but in the general well.Q Rather than work solely with the densi- 

ties, we can introduce auxiliary variables corresponding to the local potentials, so that the 
functional depends also on meson fields. The resulting DFT formulation takes the form of a 
Hartree calculation, but correlation effects can be included, if the proper density functional 
can be found. Our procedure is analogous to the well-known Kohn-Sham |35] approach in 



DFT, with the local meson fields playing the role of Kohn-Sham potentials; by introducing 
nonlinear couplings between these fields, we can implicitly include density dependence in 
the single-particle potentials. 

Moreover, by introducing the meson fields, we can incorporate the ideas of EFT. The 
exact energy functional has kinetic energy and Hartree parts (which are combined in the 
relativistic formulation) plus an "exchange-correlation" functional, which is a nonlocal, non- 
analytic functional of the densities that contains all the other many-body and relativistic 
effects. We do not try to construct the latter functional explicitly from a lagrangian (which 
would be equivalent to solving the full many-body problem), nor do we attempt here to 
construct an explicit functional using standard many-body techniques ||30|| . Rather, we ap- 
proximate the functional using an expansion in classical meson fields and their derivatives. 
The parameters introduced in the expansion can be fit to experiment, and if we have a 



1 Note that the Dirac eigenvalues do not correspond precisely to physical energy levels in the 



general case [28 
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systematic way to truncate the expansion, the framework is predictive. Thus a conventional 
mean-field energy functional fit directly to nuclear properties, if allowed to be sufficiently 
general, will automatically incorporate effects beyond the Hartree approximation, such as 
those due to short-range correlations. 

We rely on the special characteristics of nuclear ground states in a relativistic formulation, 
namely, that the mean scalar and vector potentials $ and W are large on nuclear energy 
scales but are small compared to the nucleon mass M and vary slowly in finite nuclei . 
This implies that the ratios $/M and W/M and the gradients |V$|/M 2 and \VW\/M 2 
are useful expansion parameters. Moreover, as is illustrated in Dirac-Brueckner-Hartree- 
Fock (DBHF) calculations [|37]-|39||, the scalar and vector potentials (or self-energies) are 
nearly state independent and are dominated by the Hartree contributions. Thus the Hartree 
contributions to the energy functional should dominate, and an expansion of the exchange- 
correlation functional in terms of mean fields should be a reasonable approximation. This 
"Hartree dominance" also implies that it should be a good approximation to associate the 
single-particle Dirac eigenvalues with the observed nuclear energy levels, at least for states 
near the Fermi surface ||28| . Of course, the mean-field expansion cannot accommodate all 
of the nonlocal and nonanalytic aspects of the exchange- correlation functional; however, we 
reserve a complete discussion of the limits of this approach for future work |40[ |. 

The mean-field energy functional is specified by the values of various coupling constants 
and masses. These parameters can be constrained in two ways. First, the energy functional 
provides a framework for relating model parameters to appropriate nuclear observables. This 
translates empirical nuclear properties into conditions on acceptable parametrizations. In 
practice, nuclei provide rather stringent constraints. On the other hand, the realization of 
QCD symmetries (such as chiral symmetry) in a candidate model lagrangian may impose 
relationships between parameters in the corresponding Hartree approximation to the energy 
functional. Our goal in this work is to test whether the two sets of constraints are compatible. 

We begin by analyzing a general energy functional to determine the characteristics that 
generate successful phenomenology. An essential element is a plausible truncation scheme, 
so that the functional does not entail an unrestricted number of parameters. As indicated 
earlier, our basic expansion parameters are $/M and W/M, where $ and W correspond to 
the potentials in the single-particle Dirac equations [see Eq. (§)]. We define success through 
a set of nuclear observables that should be quantitatively reproduced in any useful descrip- 
tion of ground-state properties. These observables include: 1) nuclear shape properties, such 
as charge radii and charge densities, 2) nuclear binding-energy systematics, and 3) single- 
particle properties such as level spacings and orderings, which reflect spin-orbit splittings 
and shell structure. If one recalls that the Kohn-Sham approach is formulated to reproduce 
precisely the ground-state density, and that the Hartree contributions are expected to dom- 
inate the Dirac single-particle potentials, these observables are precisely the ones for which 
meaningful comparisons with experiment should be possible. Moreover, experience has shown 
that these observables can be replaced by a set of nuclear matter properties plus constraints 
on the meson masses |I5|- PTP3"| , pi| . While the parameters of a general functional are not 
completely determined by requiring these properties to be reproduced, regions in parameter 
space that generate realistic nuclei can clearly be identified. 

The second part of our analysis is to examine specific chiral models that adopt a par- 
ticular realization of chiral symmetry at the lagrangian level and a particular mechanism 
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for spontaneous symmetry breaking. These two features restrict the accessible regions of 
parameter space, and we can simply compare with the regions determined in the first part 
of our analysis to see if there is an overlap. We stress that in this part of the analysis, we 
actually start with an explicit lagrangian and construct the energy functional by working 
strictly in the Hartree approximation; thus, the effects of correlations in these chiral mod- 
els, for example, are beyond our analysis. Nevertheless, the Hartree energy functionals are 
simply special cases of the general functional considered above, so we can apply our earlier 
results. 

The most important ingredient is a procedure that can systematically connect the nuclear 
observables and the model parameters. In previous analyses, the mean-field parameters have 
been determined either by fitting nuclear matter equilibrium properties plus some input 
from finite nuclei (e.g., the rms radius of 40 Ca) or by optimizing the predictions of selected 
experimental observables across the Periodic Table (using a \ 2 measure, for example). These 
procedures are useful if the goal is to find a good description of nuclei for a given model, but 
not if we want to rule out a class of models. In particular, it is difficult to search through 
the parameter space to determine where the boundaries of acceptable parametrizations lie. 
This severely limits qualitative insight and solid conclusions. 

The first approach was followed in our previous studies of linear sigma models: we 
explored the chiral-model parameter space by fitting to nuclear matter equilibrium properties 
and then by calculating finite nuclei for what we hoped was an exhaustive range of parameter 
sets. In the present work, we invert this analysis by adapting the approach developed by 
Bodmer ]18| , which enables us to use basic nuclear properties, derived from the observables, 
to solve directly for the model parameters. Since there are only a few basic properties that 
provide meaningful constraints, the resulting conditions are undercomplete for the most 
general functionals, and there are families of acceptable solutions defined by regions in the 
parameter space. Specific models, however, introduce relations among some parameters and 
set others to zero, so that in the accessible region of parameter space there may be no 
acceptable descriptions of nuclei. Thus one way a model can fail is if there is no overlap 
between the accessible and desirable regions of parameter space. 

There is another way a model can fail, which arises from the nonlinear meson interac- 
tions in the energy functional and in the resulting meson field equations. For some models 
and some parameter sets, one can find a nuclear ground state that has the appropriate 
properties at equilibrium, but this ground state exists only over a narrow range of nuclear 
densities. If the ground state fails to exist at densities that have clearly been obtained ex- 
perimentally (say, slightly higher than normal nuclear matter density), the disappearance of 
the model ground state must also be considered a failure of the model. We will make these 
considerations more precise as we proceed. 

The results of the present analysis solidify our previous conclusion |l7j that chiral 
hadronic models built upon the conventional linear sigma model cannot reproduce observed 
properties of finite nuclei in the Hartree approximation. The main problems in these linear 
sigma models arise from the nonlinear terms in the "Mexican-hat" potential, which serves 
to precipitate spontaneous symmetry breaking, and from the large pion coupling. The con- 
sequences are a scalar meson mass that is too large and density-dependent forces (in the 
form of scalar self-couplings) with the wrong systematics. These problems are not remedied 
by including one-baryon-loop vacuum corrections or by adding parameters that generalize 
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the model. 

On the other hand, chiral models that are phenomenologically successful can be con- 



structed. As an example, we discuss a model that was introduced in a recent paper J5B 
which features a nonlinear realization of chiral symmetry. 

For a model to be viable, it should also exhibit "naturalness" in the fitted parameters. 
Here, naturalness implies that the coefficients of the various terms in the energy functional, 
when expressed in appropriate dimensionless form, should all be of order unity. This means 
that we can anticipate the approximate magnitude of contributions to the energy functional 
(at least up to moderate nuclear densities) and thereby motivate a suitable expansion and 
truncation scheme for the functional; if the coefficients are natural, the omitted terms will 
be numerically unimportant. The naturalness of our effective field theory also implies that 
one should include all possible terms (that is, those allowed by the symmetries) through a 
given order of truncation; it is unnatural for some coefficients to vanish without an appro- 
priate symmetry argument. One consequence is that nonlinear interactions of the vector 
field should play an important role in providing a good and natural fit to nuclei, and we 
indeed find that this is the case. We stress this point because nearly all investigations of 
relativistic mean-field phenomenology to date pT|pT| - p^j26| , [27[] have arbitrarily restricted 
consideration to scalar self-interactions only.Q The quartic vector self-coupling has been 
included in phenomenological calculations only recently [ TEj, E , ^| . 

There are many studies in the literature that have characterized the nature of phe- 



nomenologically successful mean-field models, including Refs. [21-24,26,27]. We echo many 



of the common conclusions, such as the necessity for a small nucleon effective mass 
(0.58 ^ M* I'M < 0.64) at equilibrium density and for a scalar meson mass of roughly 
500 MeV. Nevertheless, since previous analyses typically determined the model parameters 
through a fitting procedure that simply calculates nuclei repeatedly until a "good fit" is ob- 
tained, connections between nuclear observables and the resulting parameters are somewhat 
obscure. Through our analysis, we try to demystify the connections between the observables 
and the parameters, so that the structure of the resulting energy functionals can be better 
understood. 

The original motivation for this work was to expand upon our previous analyses of linear 
chiral models |17| and thereby solidify our earlier conclusions. In carrying out this program, 
however, we were forced to extend our mean-field machinery in several different directions. It 
is important to emphasize these new aspects, which are more general than the chiral-model 
analysis, and which in fact comprise the most important results in this paper. First, by start- 
ing from an energy functional containing baryons and classical scalar and vector fields, we 
have a framework that goes beyond the simple Hartree approximation. This follows because 
we can interpret the analysis in the context of density functional theory and an associated 
(relativistic) Kohn-Sham approach. Second, although it is possible to vary the parameters 



2 Some previous analyses were based on renormalizable scalar-vector theories, for which such 
restrictions are appropriate. If renormalizability is abandoned, however, as in an effective field 
theory, there is no reason to omit vector-vector and scalar-vector interactions. (See the comments 
below on vector self-interactions, spin content, and causality.) 
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in the energy functional until desirable properties of nuclear matter and finite nuclei are 
obtained, it is more efficient to invert the field equations and express the model parameters 
directly in terms of the desired observables. This allows for a systematic investigation of the 
parameter space and clarifies the relationships between observables and model parameters. 
Finally, by applying ideas from effective field theory, such as the importance of naturalness 
and a suitable expansion scheme, to the construction of the ground-state energy functional, 
we find that it is important and necessary to include all allowable terms in the functional 
at the chosen level of truncation. This modern view of relativistic quantum field theory f34 ] 



generalizes earlier approaches based on renormalizable theories H3M3J20J2^] and allows for 



a more unified discussion of relativistic approaches to the nuclear many-body problem. 

The paper is organized as follows: In Section |J, we consider a general energy func- 
tional (not manifestly chiral) and characterize phenomenologically successful models. In 
Section |T|, we show how empirical properties of nuclear matter and finite nuclei can be 
used directly to determine or to constrain the parametrizations of candidate energy func- 
tionals. We also illustrate how to map out the parameter space. In Section [TV|, we specialize 
the analysis to models built on the linear sigma model with a "Mexican-hat" potential and 
demonstrate their generic failure at the Hartree level. In Section |V|, we consider a chiral 
model that successfully describes finite nuclei. Section [VT] has some additional discussion, in- 
cluding the general strategy that justifies the expansion and truncation scheme. Section |VTl 
is a summary. 

II. ENERGY FUNCTION ALS AND NUCLEAR PROPERTIES 

Our first goal is to construct an energy functional whose extremization accurately repro- 
duces bulk and single-particle nuclear observables. Experience with relativistic mean-field 
models shows that scalar and vector mesons with appropriate masses lead to NN interac- 
tions with the desired ranges, and that by introducing nonlinearities in the meson fields, 



we can include density dependence in these interactions pl|p!9| , p!8| . We assume that the 
nuclear ground states have good parity and are invariant under time reversal, which implies 
that a classical pion field does not appear^ plj. Furthermore, we retain only the valence 
nucleons explicitly. This does not mean that we are neglecting pionic or vacuum effects, but 
rather that they are implicitly contained in the coefficients of the energy functional. This 
is discussed for one-loop vacuum contributions in Ref. [36| and in Sec. 0; a more general 
discussion will be given elsewhere [HD . 



As emphasized in the Introduction, while our energy functional can be interpreted as 
arising from a hadronic lagrangian treated in the Hartree approximation, it can also be 
interpreted as an approximation to a general density functional that incorporates all corre- 
lation effects. This approach is particularly compelling because of the dominance of Hartree 
effects in the relativistic approach. This claim is supported by Dirac-Brueckner-Hartree- 
Fock (DBHF) calculations, which indicate that exchange terms and short-range correlations 



3 When nuclei are studied using mean-field models, time-reversal invariance is broken ex- 

plicitly, which leads to pion mean fields. 
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do not significantly change the size of the Hartree self-energies nor introduce strong state 



dependence (at least for occupied states) |^6| , |37| . By determining the model parameters 



from finite-density bulk and single-particle observables instead of from NN scattering, and 
by explicitly allowing meson nonlinearities that generate additional density dependence, we 
automatically include the most important effects of correlations. In fact, one can choose 
the nonlinear meson parameters so that the mean fields reproduce the scalar and vector 
self-energies obtained in a DBHF calculation [ 12 ]. 



We start with a mean-field energy functional for spherically or axially symmetric ground 
states that generalizes the functional used in Ref. P^| . It depends on a set of single-particle 
Dirac wave functions U u (x) (labeled by quantum numbers v) for the occupied valence orbitals 
and on classical, static meson fields 0(x) and Vo(x). It can be written as 

„ / occ 

E[{U„}, <P, V ] = J d 3 x C/t( x ){-ia-V + ^ M ~ + 0v^o(x)R(x) 

+ i{[V0(x)] 2 + m b 2 [0(x)] 2 ]} - i{[W (x)] 2 + m 2 [V (x)] 2 } 
+ !«#(x)] 3 + ^A[0(x)] 4 - ^C[^Vo(x)] 4 

+ a-[g v V (x)Mx) - ia'[^y o (x)] 2 [0(x)] 2 ) , (1) 



subject to the constraint 



d 3 x^(x)^(x) = l, (2) 



for all occupied states. Observe that the functional has been truncated at quartic terms in 
the fields and at quadratic terms in gradients of the fields; this truncation will be justified 
below. 

To realistically describe finite nuclei, the functional must be extended to include rho 
mesons, pions, and photons. In the present discussion, however, we focus on nuclear ground- 
state properties that primarily constrain isoscalar physics. Thus the Coulomb contribution 
and a correct reproduction of the bulk symmetry energy (« 35MeV) are sufficient for our 
purposes, and this requires only the conventional extensions to include neutral rho mesons 
(denoted by bo) and the Coulomb field (see also Ref. |p0| ): 

E[{U U }, 0, V , bo, Ao] = E[{U U }, 0, V ] + 

, OCC 1 1 

Jd 3 x (J2 ^(x){-^r 3 6 (x) + e-(l + r 3 )A (x)}^(x) 

" ^{[V&o(x)] 2 + m 2 [& (x)] 2 } - i[VA (x)] 2 ) . (3) 
For simplicity, we suppress these isovector terms in the sequel, although they are included 



in the numerical calculations of finite nuclei. As noted above, the chiral models of Sees. [TV 
and [V] have no explicit pion field contribution in this formulation. 

The meson fields are determined by requiring the functional to be stationary with respect 
to their variations; this yields (after partial integrations) the meson field equations 
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(V 2 -m s 2 )0(x) = -g s J2U u (x)Uv(x) 

V 

+ i*#(x)] 2 + ^A[0(x)] 3 + (a^ - a>(x))[^y (x)] 2 , (4) 

occ 

(V 2 -m 2 )V (x) = -g v Y,Ul(x)U v (x) 



V 



+ W v lVo(x)T - 2«-^ 2 \/ o (x)0(x) + a'g 2 v V (x)[<j>(x)] 2 . (5) 

The single-particle orbitals are determined similarly. The constraint equation (Q) is imposed 
for each orbital U u by using a Lagrange multiplier e u , which we identify as the energy 
eigenvalue of the Dirac equation: 

{-za-V + (3[M - <? s 0(x)] + g w V (x)}U v (x) = e v U v {x) . (6) 

The energy is minimized by solving these equations self-consistently for the N and Z lowest 
eigenvalues to determine the occupied neutron and proton states. 

This energy functional, which provides an approximation to a general density functional, 
could also be derived as the one-baryon-loop energy from the lagrangian density 

C = - 9,1^ - (M — g s <f>)}il> + 1 -{8^<P - m s 2 2 ) 

- \{d,V v - d v V,f + l -mlV,V» - ^ - 1A0 4 
+ ^UiV^f - agl^Vy + ^a'gXvy 2 , (7) 



by following the discussion in Ref. fS^j.f] Alternatively, one could construct the Hartree 
Hamiltonian from Eq. (^|) and take its expectation value in a state specified by static meson 
fields and a set of occupied single-nucleon orbitals (labeled by quantum numbers v). In this 
lagrangian formulation, it appears that vacuum contributions are neglected, because the sum 
in E runs over valence nucleons only. However, the vacuum effects can be precisely absorbed 
into the coefficients of the interaction terms, provided we keep nonlinearities to all orders 
f36fl . If we assume that a truncation at some low order of the fields and their derivatives is 
sufficiently accurate at nuclear densities (which we justify on the basis of naturalness), then 
the vacuum effects are already included with sufficient accuracy in the functional of Eq. ([[]). 

The notation and structure of the lagrangian (|7]) has been chosen to conform to previous 
usage in the literature. It includes as special cases most of the lagrangians used in mean- 
field studies. The Walecka model retains the terms with nucleon fields and just the kinetic 



4 Note that this derivation of the functional starts from a Lorentz-invariant action S = J d xC 
If we consider the energy functional as an effective functional, however, we presently know of no 
reason to exclude terms that explicitly contain the medium four- velocity u^, such as u IJ, V fj y u V u . 
This issue will be considered in a later publication. 
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TABLE I. Relationship between parameters of Eq. ([[J) or Eq. (Q) and specific models. Chiral 
models with a linear representation should have g s = g n at the level considered here; however, we 
will allow g s = g n / 1 g& with g n ~ 13.5 and gA ~ 1-26, so that the Goldberger-Treiman relation is 
satisfied. We have also taken = in writing the entries. Note that the rfs in the last two rows 
are unrelated. 



Model 


Ref. 


9v 


K 


A 


c 


a 


a' 



Walecka 
Nonlinear scalar 

Bodmer 

Chiral auj 

Boguta 

General chiral 

Nonlinear chiral 



m 



m 
m 



i 







9w 
9v 

9v 
9v 

M 
9v 





-2b 

2a 





6c 
6b 



3g s m 2 


3g 2 m 2 


M 


M 2 


3g s m 2 


3g 2 s m 2 s 


M 


M 2 


3g s m 2 


3g 2 m 2 


M 


M 2 


(3d - 8)m 2 s 


(lid 2 -A8d + A8)m 2 


dS 


(dSoY 






6m 2 




C 
C 









1 

>f 

~V9s m v 
2S^M^ 









1 

>f 




and mass terms for the scalar and vector fields. 



introduced long ago by Schiff |47| and later studied by Boguta and Bodmer 



Nonlinear terms in the scalar field were 

Although 



widely used, a long-standing worry about such models was that the quartic scalar coupling 



A is negative in good fits to finite nuclei |21|J23|J24| | 

This concern is eliminated when quartic vector self-interactions are included in the la- 
grangian or energy functional. Such a term was first introduced by Bodmer and Price 

13, IS], motivated by a search for softer high-density equations of state. Gmuca employed 
the quartic vector term to account for density dependence in the vector self-energy, as im- 
plied by DBHF calculations [12 . Previous studies of lagrangians with neutral vector mesons 
have suggested that such terms can lead to field equations with causality- violating solutions 

HQ or to the mixing of spin-one and spin- zero representations |5l| . However, if ( > 



17 



and the other vector parameters are natural, then there are no problems with causality. The 
second point is not an issue in our effective-field-theory approach; we are not describing 
elementary vector mesons, and we certainly expect mixing effects at finite density [55 . 

Equation ([?]) also includes as special cases all of the chiral models considered in Ref 
and the chiral model of Ref. |36[.[] The relationships between parameters in Eq. (|7|) 
models from the literature are summarized in Table |. 

While our immediate intent is to work with a functional that includes the chiral models 
as special cases, the proposed functional (|l|) is actually more general than it appears at first. 
One can imagine several classes of additional terms in Eq. (0): 



and 



5 The model of Ref. |3!| actually contains a logarithmic potential for the scalar field, which means 
an infinite polynomial in <j). However, for reasons discussed below, specifying the cubic and quartic 
terms determines the energy functional at ordinary densities at a level sufficient for our purposes. 
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1. Higher-order self-couplings and derivatives involving meson fields alone [e.g., 5 or 
(^%) 3 or {d^d^) 2 ]; 

2. Contact terms involving nucleon fields beyond bilinear order [e.g., (ipijj) 2 ] ; 

3. More complicated meson-nucleon couplings [e.g., g a ((f>)ipip<j)]. 

Let us consider each in turn. 

For applications of the energy functional to ordinary nuclei, contributions to the energy 
from higher-order polynomials (beyond quartic) or higher-order gradient terms (beyond 
quadratic) are numerically small, unless the coefficients are "unnaturally" large. The defini- 
tion of naturalness proposed here (and discussed further in Sec. |VTD is simple: when written 
in appropriate dimensionless form, the coefficients of all mesonic terms in E are of order 
unity. Thus we can organize our truncation by counting powers of the local scalar and vec- 
tor potentials (or self-energies) that appear in the single-particle Dirac equations, divided 
by the nucleon mass, and by counting gradients of the potentials divided by the square of 
the mass. We take as a basic principle that our functional implicitly includes higher-order 
terms, but with natural coefficients, which can be ignored in practice at ordinary nuclear 
densities and below; their small contributions can be absorbed into slight adjustments of the 
other coefficients. Thus we have the possibility of a useful expansion and truncation scheme 
in this framework. We must establish, however, after determining the parameters, that the 
highest-order terms retained do not dominate the energy, and that adding additional terms 
produces only small changes in the parameters. These considerations and the consequences 
of this approach for extrapolations to high density are discussed in Sec. [Vl|. 

Contact terms arise if we "integrate out" the scalar and vector fields (by using the 
field equations, for example). Conversely, contact terms included originally in the energy 
functional (excluding certain terms with derivatives) can be eliminated in favor of the scalar 
and vector fields, if we allow products of fields to all orders. The issue then really becomes 
one of efficiency, since in either framework one will have to truncate in practice. Based 
on the economical successes of relativistic meson-exchange phenomenology and to connect 
to specific models, we use the field expansion. We again rely on the naturalness of the 
coefficients to motivate a truncation at fourth order in the fields.^ Further discussion and 
development of these ideas will be given in a future paper J4tj . 



Models have been proposed that replace M* = M — g s (p in the lagrangian by a more gen- 



eral function M*(0) [or equivalently, that make the replacement g s — ► g a (<f>)] |5^H54|- This 
function must have a Taylor expansion around vanishing <p, but is otherwise quite general 
in principle. More generally, the Yukawa couplings of the scalar and vector to the nucleon 
fields can be replaced by polynomials in <fi and (in Lorentz-invariant combinations). The 
motivation for such models is often that these terms are needed to reflect the compositeness 
of the nucleon. However, if we allow arbitrary polynomials in the meson fields in our func- 
tional, we can simply redefine the fields to eliminate non- Yukawa couplings to the nucleon. 



Strictly speaking, in this approach m v should be a fit parameter rather than fixed at the physical 
u meson mass. We rely on resonance dominance and the insensitivity of the functional to the 
precise value of this mass to justify fixing it. 
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(For example, we can replace <ft <p where M*((f>) = M — g' s 4>.) The new functional has the 
same form as Eq. (Jj), including all orders in the fields alone. The observables are unchanged 
by the transformation, and the only question is again one of efficiency, that is, whether a 
truncation to fourth order in the fields is an accurate numerical approximation. 

Thus, by working with Eq. (P, we accommodate a wide class of effective models. To start 
our analysis, we identify the nuclear observables that we wish our functional to reproduce. 
The properties of finite nuclei fall into three major categories: 



1. Nuclear shapes. This includes the basic manifestations of nuclear saturation pREH]: a 



flat interior, a surface thickness independent of the baryon number B, and a systematic 
increase of the nuclear radius with B 1 ^ 3 . These features are seen experimentally in 
the charge radius and the charge density; the latter is most clearly analyzed for this 
purpose in momentum space (i.e., the form factor). 

2. Binding-energy systematics. The binding energies of (closed-shell) nuclei across the 
periodic table must reflect the liquid-drop systematics of the basic semi-empirical mass 
formula. At a coarse level, this involves the sensitive interplay of the bulk binding 
energy (ai), the surface energy (02), the Coulomb energy (as), and the bulk symmetry 
energy (0,4) (with a surface correction included) |55|-^]:[] 

Z 2 (N — Z) 2 

E/B - M = -a 1 + a 2 B~ 1/3 + a 3 —— + a A ^r-± . (8) 

1 6 b 4 / 3 B 2 (l + 3.28/B 1 / 3 ) KJ 

3. Single-particle properties. The single-particle potential is reflected in the ordering and 
spacing of single-particle levels, which are reasonably well identified experimentally 
(up to rearrangement effects that are comparatively small). Speaking in nonrelativistic 
language (since the Dirac equation can always be recast in "Schrodinger-equivalent" 
form), the level ordering in I implies a central potential that interpolates between 
a harmonic oscillator and square well. In addition, the spin-orbit potential must be 
strong enough to ensure the correct shell closures, but not too strong, or else the 
systematics of nuclear deformations will not be reproduced |[22|| . Finally, there is the 
energy dependence of the optical potential at low energies, which is related to the 



vector self-energy at the relativistic mean-field level [20 



One could add other properties to this list. In many cases, however, these additional prop- 
erties are strongly correlated with the features listed above and so are reproduced without 
imposing further conditions. 

Using properties of nuclei as constraints, however, is difficult in the sort of direct analysis 
we seek. We will therefore extrapolate from the systematics of actual observables (such as 
binding energies, rms radii, and spin-orbit splittings) and impose conditions on the func- 
tional in nuclear matter near the equilibrium density. First, nuclear matter must exhibit a 



There are much more sophisticated mass formulas incorporating additional physics, but this is 
the appropriate level of detail for our purpose. 
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particular equilibrium density (po) and binding energy (eo) within a fairly narrow range. 
These conditions are supported by calculations in relativistic models that are fit directly to 
properties of finite nuclei, which consistently predict similar p an d eo when extrapolated to 



infinite nuclear matter 23.26 



We can also consider derivatives of the energy with respect to the density evaluated at 
equilibrium. The second derivative or curvature determines the compression modulus (Kq). 
Even though the range in Kq consistent with nuclear properties is fairly broad (approxi- 
mately 200-300 MeV), this provides significant constraints on mean- field models. The most 
direct effect of the compressibility on ground-state properties is through the surface energy 
and is therefore manifested in the energy systematics. The surface energy is also correlated 
with the scalar mass. 

The importance of higher derivatives of the energy is unclear at present. Recent work has 
suggested that the ratio of the "skewness" (related to the third derivative) to the compression 
modulus is well determined by nuclear monopole vibrations [^,0 • We do not use skewness 
to constrain the energy functional in this investigation, but we comment on the correlation 
between skewness and ground-state properties in Sec. |VT[ Note that one can easily extend 
the analysis described in Sec. |T| to exploit additional constraints of this type. 

To correctly reproduce empirical single-particle energies, the most important ingredient 
beyond the depth and shape of the effective central potential is the strength of the spin-orbit 
interaction. In relativistic models, the spin-orbit splittings are highly correlated with the 
nucleon effective mass M*. An upper bound on M* is given by the reproduction of the 
shell closures in heavy nuclei and the splittings of spin-orbit pairs in light nuclei, while the 



reproduction of deformation systematics gives a lower bound [2^]. In conjunction with the 
other constraints!] Mq/M must be in a fairly narrow interval around 0.6. 

The symmetry energy is determined by the observed energy systematics [through Eq. (]8|)], 
the relative energies of proton and neutron single-particle levels in heavy nuclei |19| , and iso- 
tope shifts in neutron-rich nuclei P6"| , |6"T| . These imply that ~ 35 MeV, with an uncertainty 
of several MeV. 

Finally, there are important restrictions on the ranges of the attractive and repulsive 
interactions, or equivalently, on the momentum dependence of the scalar and vector fields. 
These conditions do not impact nuclear matter (since the momentum transfer is zero at 
the mean-field level) but are critical in finite nuclei. The major constraint is on the scalar 
mass, which is consistently found from direct fits to be in the vicinity of 500 MeV (when 
the vector mass is taken to be the physical u mass). Significantly larger scalar masses 
produce unobserved oscillations in the charge density and too steep a surface in the effective 
central potential, which leads to incorrect level orderings. We will assess the importance 



We will use the convention that properties and fields evaluated at equilibrium density are denoted 
with a subscript "0" . 

9 Spin-orbit strength is not the only source of constraints on M*. Reinhard et al. concluded that 
Mq/M must be roughly 0.6 from detailed fits to finite nuclei that did not include the spin-orbit 
splittings as input plf . 
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of this constraint for chiral models (and the success of the nuclear matter constraints) by 
calculating finite nuclei after the nuclear matter analysis. 

One may ask: at what level of accuracy should we require nuclear observables to be 
reproduced? The answer depends on the planned application of the model and the goals of 
the analysis. For comparison with nuclear scattering experiments, accurate nuclear densities 
and wave functions are required, but total binding energies are not particularly important; 
in contrast, for studying fission barriers or the properties of nuclei near the "drip lines", 
energy systematics are paramount. Since our goal is to exclude some classes of mean-field 
models, we will impose less stringent constraints. The nuclear matter equilibrium density 
Po and binding energy cq are rather well known, and small variations are not important to 
our conclusions. We therefore fix 

p = Pb0 = 0.1484 fm~ 3 (Jfejj = 1.30 fin -1 ) , (9) 
g 

e = — - M = -16.1 MeV = -a x . (10) 
Po 

We choose the ranges of acceptable K$ and Mq based on models that have been directly 
optimized for finite-nucleus observables, allowing considerable uncertainties. Specifically, we 
consider 

180 MeV < K < 360 MeV , (11) 



0.58 < M*/M < 0.64 . 



(12) 



If one focuses on a subset of the nuclear observables, one can stray outside these bounds 
and still obtain good results. For example, the original Walecka model generates quite 
acceptable nuclear densities despite its high compression modulus, which yields unfavorable 
energy systematics |19 



III. CONSTRAINING NUCLEAR MATTER ENERGY FUNCTION ALS 

Given the conditions on the energy functional for nuclear matter [namely, eo, Po, Kq, and 
Mq in Eqs. ©-(0)]; our g° a l is to solve for the parameters of the functional as functions of 
these input "observables." In general, we can determine four of the parameters from these 
four inputs, and the additional coefficients will then specify the space of acceptable models. 
To rule out a class of models, as we will do in Sec. [TV], we must show that it is impossible 
to find a set of coefficients that fall within the bounds of acceptability. 

In uniform nuclear matter, the energy functional of Eq. (JTJ) simplifies to the energy 
density function 

£{4, V ; p B ) = g v V p B - \m\V 2 + g 2 v V 2 (a-<t> - l -a>cf) - ±(gX 

+ \nV + + + 7£jl / V P <P 2 + M * 2 ) 1/2 • < 13 > 
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where 



M* = M-g s 4>, (14) 
and the Fermi momentum is defined through 

The spin-isospin degeneracy 7 is 4 for symmetric nuclear matter. The coefficients k and A 
may be constrained by chiral symmetry (as in the linear sigma model), but are free at this 
point. Note also, that if we choose to work with a renormalizable lagrangian and include the 
one-baryon-loop vacuum energy A£ vac (either chirally or not), this is equivalent to redefining 
k and A as well as the coefficients of higher powers of <ft (which are not shown). As before, 
the parameters and their special values in particular models are given in Table |I[ 

It is useful to change variables in the energy density. First we switch from <ft and Vq to 
the mean-field nucleon self-energies (or optical potentials) $ and W (following the notation 
of Ref. 0): 

$ = g B <f> , W = g Y V . (16) 

The meson field equations are now determined by extremizing £ with respect to $ and W . 
As discussed earlier, the natural expansion parameters for the finite-density functional are 
W/M and $/M. The quadratic terms in Eq. ([13]) imply that the couplings g s and g v do not 
appear individually in nuclear matter, so we rewrite the energy density in terms of variables 
that absorb the couplings. We introduce the ratios of meson couplings to masses 



A = lb. 



(17) 



together with the scalar self-couplings 



77-4, A=4> (18) 



S S 3 ' " Ss 4 



and the scalar-vector couplings 



_ _ a _, _ a! 

9; 9; 



(19) 



The correspondences between conventional model parameters and the newly defined param- 
eters are summarized in Table |I[ The nucleon mass M = 939 MeV and the vector meson 
mass m v = m w = 783 MeV will be fixed at their experimental values. 
In terms of the new variables, the energy density can be written as 

W- p B ) = £ v$ ($, W; p B ) + U($) + £ fe ($; p B ) , (20) 

where 
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TABLE II. Relationship between parameters of Eq. (|T^) and specific models. The quadratic 
coupling is c s = g s /m s . Chiral models with a linear representation should have g s = at the level 
considered here; however, we will allow g s = g w /gA with g n « 13.5 and gA ~ 1.26, so that the 
Goldberger-Treiman relation is satisfied. We have also taken m n = in writing the entries. Recall 
that the 77's in the last two rows are unrelated. 



Model Ref. 


Cy 


K 


A 


C 


a 


a 1 


Walecka Q 


m v 


















Nonlinear scalar @] -2| 6c Q Q 



9* 



Bodmer © £ | » §^ 

Chiral.. © £ _^ V 



General chiral @ £- ( ^ ^ 

Nonlinear chiral M {M ~ 8) .. (j^ligj+ig) /■ z!L n 

Nonlinear chiral ^ mv ^^-^ (d&So)^ C 2{g s S»)c*M 



£ v$ ($, VT; p B ) = l^p B - ^I^ 2 + W 2 (aM§ - l -a'&) - 



(21) 



V ; 2c s 2 6 24 



(22) 



and £fc is the "kinetic energy" (including the nucleon rest mass): 
£ k (^p B ) = | fcF d 3 Pv /p 2 + M* 2 

7 



8tt 2 



A; fj E; 3 - —M* 2 k F E* - -M*Hn(^^ 
2 F 2 V M* 



with 



e; 



kl + M* 2 



(23) 



(24) 



This division is useful because it isolates the W dependence and the parameters in U (c 2 , 
7c, and A), for which we will solve. We could also solve for a different set, as in Sec. |TV]. 

The energy per particle (S / p B — M) and pressure p are determined by the seven free 
parameters c v , c s , k, A, a, a', and (. Any two models with the same values of these 
constants will give identical results for nuclear matter, so these are the relevant parameters 
for testing our constraints, rather than those in the original energy functional. 

Let us summarize the conditions on £ as a function of W, and p B (T^|. By evaluating 
these conditions at equilibrium, we will be able to solve for the parameters in terms of the 
input data e , Po, K , and Mq. We observe immediately that if Mq is given, then $ follows 
trivially from 
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$ Q = M - M * . 



(25) 



It is therefore convenient to define the equilibrium values 



(26) 



and 



us = u"(%) 



dE7($) 



d$ 



d-o c s 2 2 6 " 



d 2 C/($) 



d$ 2 



1 



+ k$ + -$ 2 



(27) 



(28) 



The conditions on S are as follows: 

1. At any density p B , £ is stationary with respect to $: 

(I) = »■ < 29 > 

Carrying out this derivative on Eq. ( ^0|) and evaluating the result at equilibrium yields 

p s0 - W*[aM - a'%] - U' Q = , (30) 



where 



with 



PsO = Ps($o,Po) 



(31) 



Ps($,Pb) 



d§ / w, PB 



7 



(2tt) 3 



d 3 P 



M* 



Vp 2 + M* 2 



M* 



(32) 



2. The vector field appears only in £ v $ and is determined by the equation of constraint 

, (33) 



/ dS_\ / ggyg 

\dw)$ lPB ~ \dw y$, PB 



which produces (at equilibrium) 



P„ - -^Wo + W [2aM$ - a'<S> 2 ] - |w o 3 = 



1 



6 



(34) 
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3. The binding energy per particle at p B = po is — eo, so 

S = (eo + M)p • 



(35) 



4. The Hugenholtz-van Hove theorem follows by calculating the chemical potential p 
from £ in Eq. (EH): 



d£ 

dp B 



V<9p E 



B / * 



W + A /fc2 + M 



*2 



(36) 



Note that we need only take (d£/dp B )^ t w to obtain this result, since £ is stationary 
with respect to $ and W from Eqs. (E3) and (|3"3"D. The mean-field energy density is 
thermodynamically consistent, so the pressure follows from the first law of thermody- 
namics (as can be explicitly verified by taking derivatives): 



p = -£ + p B p = -£ + p B W + p B E* 



(37) 



The condition for equilibrium is p = at p B ^ 0, so combining Eq. 
evaluated at equilibrium yields 



W = e + M- J(k F0 ) 2 + Mp 



with Eq. (|37|) 



(38) 



This quite general result provides a direct connection between Wo and the inputs e , 
po, and M* 



5. The compression modulus K follows from 



,d 2 (£/p E 



dp : 



Pb 



d 2 £ 



dpi P 



2-2-, 



(39) 



so that at equilibrium we have 

d 2 £ 



K = K(p ) = 9 Po 



9p 



dpi 



7T 



+ 



dW 



2k F0 E* dp 



fdW 



-fo.Po + V 9$ 



M* \ d$ 



*o,po E* J dp 



P() 



Carrying out the derivatives using the field equations^ produces 



7T 



v 



( E#- - V /( 



fO-C'fO 



^'-Pso + ^vo-a'^ 



(40) 



(41) 



10 This is most simply done by considering TV to be a function of $ and p B . 
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where 

2 



... (Pso - U' Q ) = 2W (aM - s'$ ) , (42) 
Wo 



v = 



Wp 1 

p + RWo 1/4 - (2a$ M - a>& ) + ICW* ' 



(43) 



and 



P'sO 



dp s 



dp s 



dM* 



(44) 



with 



dps 



M* vr 2 



/bp n / /up ~t~ £/ F 

— m' 



El 



M* 



(45) 



We can now systematically solve for c 2 and the three parameters in £/($) using the 
four inputs eo, po, ^o, an d Md- This leaves as free parameters (, a, and a', which label 
potentially viable solutions. We note immediately that Mq gives us <3>o from Eq. ( p5|) and 
then Wq from Eq. fl3"8p. We can now solve Eq. ( |54] ) for c 2 , since given £, a, and a ', everything 
else is known: 



Wn 



v p + W [2aM% - a%] - \CW$ ' 
Next, we express c 2 , 7c, and A in terms of C/q, Z/q, and Z7q by inverting Eqs. (|26"|)-(p8"l):pl 



(46) 



K 



A 



2 




- 3$ ^ + 6f/ 




+ 5$ t/ - 8U 




6 




- 2$ ^ + 3f/ 


24^0 



(47) 
(48) 
(49) 



The final step is to relate Z7 , £/q, and Uq to the input observables. After determining 
$ ) Wo, and cv, we can calculate £ v $ and £k at equilibrium. Then Uq = f/($ ) follows from 
Eqs. flSj) and (0): 



n At this stage, the form of U(<&) is not critical, and if U ($) has three unknown parameters (but 
not necessarily a quartic polynomial), one can just calculate Uq, Uq, and Uq and invert to find the 
parameters in terms of poj Mq, and Kq [pj[]. 
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U = £ - £ v *($o, W ; po) - £ k (%; p ) 

= (e + M)p - £ v$ ($ , W ; p ) - £ k (%; p ) . (50) 

Here e , Po, Mq, and ( are given, £fc($oj Po) follows from Eq. (^) at the given p and Mq, 
and we use Eqs. ([H]) and (|34D to get £ v $(<3>o, Wo; po), with the result 



£ v *(%, W ; po) = ~WVo + • (51) 



We can evaluate U' using Eqs. (p0|) and (|32l): 

f/o = Pso- »?[sM-a'$o] (52) 
and then solve the K equation (HT1) for £/q in terms of known quantities: 



[ - v / ( 

US = P' s0 - fro + a <Wl + ^ . (53) 

If we use all four inputs e , Po, Mq, and fT to determine the parameters c^, Cg, ~K, and 
A, we can find solutions parametrized by the remaining free couplings, which are (, a, and 
a' in the general case. It would seem that we can always find a solution. How, then, can we 
ever exclude a model? There are several possibilities. First, the model may be less general, 
so that there are not enough free parameters. For example, the Walecka model has only 
two free parameters (c^ and c^); thus, only e and p can be specified, and Mq and K are 
predicted. There is a transcendental equation for Mq to be solved in this case [see Eq. (p5| . 
below], and the solution may not fall within the bounds specified by Eq. ([12]). Another 
possibility is that the solution manifests unphysical conditions (e.g., < or c\ < 0). 

A third possibility is that the parameters allow multiple solutions to the meson field 
equations, which may correspond to pathological or abnormal solutions (e.g., a Lee- Wick 
solution [|J in the sigma model). These alternative solutions extremize the energy functional 
and may have a lower energy than the normal solution (which satisfies the phenomenological 
input conditions). As discussed in the next section, we will not necessarily exclude a model 
in this case unless the normal solution, which exists by construction at p , disappears at a 
slightly higher density. If this occurs, we must consider the parameter set unacceptable. 

To investigate the latter possibility, we must always check for multiple solutions to the 
$ and W equations [(EDI) and (133)] after determining the parameters. Equation (R¥) is easy 



to analyze, since it is a cubic in W with known coefficients (for given $). For example, one 
can show that with a = a' = 0, it is necessary that ( > to avoid abnormal nuclear matter 
solutions with lower energy than the normal solution. We illustrate what can happen with 
multiple solutions in Fig. ^ for a = a' = and ( = —0.03. The solid line is the normal 
solution, with a minimum at eo and po- The other curves arise from different solutions to 
the W and $ equations (three each!). The normal solution disappears at a density slightly 
above po, and only pathological solutions remain at higher densities (not shown). 

A final possibility is that the parameter set is unnatural, that is, the coefficients of the 
highest-order terms are too large. These terms then dominate the energy functional at 
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FIG. 1. Binding-energy curves with £ = —0.03, a = a' = 0, and all other parameters determined 
from eo, po, -Mq, and -Ko- All of the curves extremize the energy functional. The normal solution (with 
equilibrium binding energy — 16.1MeV at k F = 1.3 fen ) is the solid line. Only four of the six solutions are 
shown. 



densities near nuclear matter equilibrium, which implies that our truncation is unjustified, 
and the entire framework breaks down. 

The Bodmer method for determining the model parameters from information about 
nuclear matter equilibrium can be conveniently implemented in a spreadsheet or symbolic 
manipulation program. It can be summarized as follows: 

1. Specify e , Po, K , Mq, and also (, a, and a'. (Here we have chosen to solve for c^, 
Cg, k, and A.) 



2. $ follows from Eq. (|25|) and W from Eq. (|3 

3. Determine from Eq. (f46|). If m v is specified, then g\ follows immediately from the 
definition in Eq. (17). 

4. Compute U from Eqs. (|5|), and W Q from Eqs. © and (||); and Ug from 
Eqs. ©-© and Solve for Cg, k, and A from Eqs. (f47l) -(f49l). If the scalar mass 
m s is also specified, one can obtain g s , k, and A from Eqs. (|TT1) and (^8|). 

5. Find all solutions to the scalar and vector equations. (It is easiest to first eliminate 
W from the $ equation by directly solving the cubic W equation.) Check that the 
normal solution is acceptable within a prescribed density range. 

As a starting point for our discussion, we set a = a' = 0, so that we consider the same 
class of models studied by Bodmer [0. A useful procedure is to fix eo and po and then 
to scan the (Kq, Mq) plane, with each point producing a family of parameter sets, since 
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there are five free parameters and only four constraints. We can then establish boundaries 
corresponding to special values of the parameters. 
With a = a' = 0, the energy density reduces to 

£ = Wp B + £ fc ($; p B ) - ^W* + ^ + ^ + lA$ 4 - ±(W* . (54) 

This implies that 

• ( > lowers the energy (attractive), 

• k, A > raise the energy (repulsive), 

• k, A < lower the energy (attractive). 

If k — A = (no scalar nonlinearities) and £ = (no vector nonlinearities), then we obtain 
the original Walecka model. In equilibrium, we can set £ from Eq. (|54"|) equal to (eo + M)p , 
and eliminate Wq using Eq. ([8J), q? using Eq. (^), and using Eq. fl3~4]). This leaves a 
transcendental equation for Mq: 

(e + M) Po + V( fc Fo) 2 + Mo 2 Po — (M — M*)p s0 - 2£ k (M*; p ) = . (55) 



For the values of p an d e in Eqs. (§) and (|10D, the solution is Mq/M = 0.538, which then 
implies K = 557 MeV from Eq. (f41~D. This point is denoted as (A) in Fig. |. 

If C = but k 7^ and A ^ 0, it is possible to reproduce the desired Mq and K with 
[/($) alone, which is well known. In this case, k > and A < 0. The situation can be 
displayed graphically in the following figures. In Fig. 0, we plot the locus of points for which 
k = or A = 0, as a function of both Mq/M and K . This lets us separate the regions 
of parameter space corresponding to different signs of k and A. The boxed area shows the 
acceptable valuesof the input "observables" (0.58 < M*/M < 0.64, 180 < K < 360 MeV). 
Note that when A < 0, the energy is unbounded below for large this is not necessarily a 
problem, however, since our truncated functional is only valid for small 

If we now allow ( ^ (recall that ( must be positive), then as ( increases, we get 
additional attractive contributions to the energy, which leads us to increase A. Figure |3| 
shows lines of A = for selected values of c1(Wq/6. [This combination is the ratio of the 
cubic and linear terms in Eq. ( |3"4"D and is thus the relevant measure of vector nonlinearities.] 
For large enough (, we can achieve acceptable Mq and K with A > 0. We can see how 
this works by examining the expression for £ in Eq. (|5*4"D. Since both $o and Wo are 
approximately proportional to p B , a phenomenologically acceptable £ clearly requires an 
explicit dependence cx p B with a negative coefficient. If ( = 0, we must set A < to achieve 
this. By setting ( > 0, however, we can increase A and obtain essentially the same result. 
Nevertheless, we must be careful not to let the quartic nonlinearities get so large that they 
dominate the energy, for this would violate our assumption of naturalness. 

If we hold K = X = a = a' = and let ( > (vector nonlinearities only), then Mq/M 
decreases and K Q increases relative to the Walecka model, as indicated by the line in Fig. |j. 
Thus adding ( alone makes the resulting nuclear matter properties worse; to return to the 
desired region of Mq and K Q , we must restore the scalar nonlinearities. Furthermore, when 
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FIG. 2. M */M vs. K Q with C = a = a' = 0. The lines show the k = and A = boundaries. The box 



encloses the desirable values of Mq/M and Kq. Model B is from Ref. |17j (see Table III) 
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FIG. 3. Mq/M vs. Ko with a — a' — 0. The lines show the A = boundaries for a range of values of 
Cy(Wo/6 from 0.0 to 0.2. This combination is a measure of the vector nonlincarities. The box encloses the 
desirable values of Mq/M and K . 
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FIG. 4. Mq/M vs. Kq with a = a' = and n = A = 0. The line shows possible solutions as a function 
of (c^Wq /6. The C = point is the Walecka model. The box encloses the desirable values of Mq/M and 
K . 



one increases £, one must compensate by increasing A, so that an acceptable description 
is obtained only when both terms are included. This conclusion is consistent with our 
definition of naturalness and is also supported by the density dependence of the scalar and 
vector self-energies found in DBHF calculations 



One can further explore the space of acceptable parametrizations by letting a ^ and 
a' 7^ 0. Although the enlarged parameter space allows for a wide variety of acceptable 
solutions, this additional freedom is not needed to satisfy the phenomenological constraints 
considered here. Nevertheless, two observations can be made. First, including the a and a' 
parameters allows one to adjust the density dependence of the effective scalar and vector 
meson masses (defined by diagonalizing the matrix of appropriate second derivatives of the 
energy functional), which may provide useful constraints in the future if concrete empirical 
information becomes available. (Note that one must be careful of the distinction between 
longitudinal and transverse masses when trying to compare with experiment.) These addi- 
tional parameters may also allow for improved fits to DBHF self-energies or to self-energies 
obtained in more sophisticated nuclear matter calculations. Second, our definition of nat- 
uralness says that there is no reason to omit these scalar-vector couplings, unless there is 
some (as yet unknown) symmetry principle that forbids them. In the interest of brevity, we 
leave an exploration of the expanded parameter space to a future investigation. We turn 
instead to specific chiral models, considered as special cases of the full energy functional, 
and apply a variant of the Bodmer analysis described above. 
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IV. ANALYSIS OF LINEAR SIGMA MODELS 



We now adapt the general analysis of the previous section to consider mean-field models 
built on the linear sigma model with a light scalar meson. There is a long history of such 
models, starting with the lagrangian considered by Kerman and Miller [p|,|20|: 

C = #7 M <9 M - SVT^V - cjw(<7 + ry 5 T-7r)]V> 

+ \{d^a + d^-d^) - -\{a 2 + 7T 2 - v 2 ) 2 

- \{d^V v - d u V^ 2 + \mlV^ + ea . (56) 

Because of the "Mexican-hat" potential, v 7^ implies that a acquires a nonzero vacuum 
expectation value, which generates masses for the nucleon and scalar meson (and pion). This 
also fixes the self-couplings of the scalar field <fi that is the deviation of a from its vacuum 
value. Here we study descriptions of nuclei that arise from generalizations of this model. We 
stress that we are now exploring the phenomenological consequences of definite lagrangians 
treated in the Hartree approximation. We are able to use the effective density functional 
analysis developed in the previous sections because the form of the energy functional is the 
same. 

The general extension of the linear sigma model (|56[) considered in Ref. [17]] is based on 



the lagrangian (after shifting the scalar field a from its vacuum value to define 0) 
C = #7^ - g^V - (M - gA) - ig^r-ir]^ 
+ i(0„ W - m s V) + \{d^-d^ - mW) 

+ - mlm 2 + n 2 ) - ^(m 2 - mlM + n 2 ) 2 

- \{d,v v - d v v,) 2 + \mlv^ + ^Cgtiv^f 

, MVpVt + lv 2 g%v^4> 2 + * 2 ) . (57) 
g-K j 

Here we have eliminated A, v, and e in favor of M, m s , and m n . We will take m n = 
in the sequel.0 This is the most general lagrangian with nonderivative couplings through 
dimension four that is consistent with linear chiral symmetry (for m v = 0). 




When solving chiral mean-field models, one must be aware of the possibility of bound "anoma- 
lous" solutions, in which the scalar field interpolates between the minima of the effective potential. 
Including a finite pion mass increases the energy of the anomalous solution so that it is above 
the energy of the normal solution. Hence we can simply ignore these anomalous solutions and set 
m n = 0. (The pion mass has negligible effect on "normal" solutions.) 
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Note that chiral symmetry in the linear realization implies that the scalar-nucleon 
Yukawa coupling constant g s is equal to the pion-nucleon coupling constant g w , which is 
known experimentally. Moreover, the nucleon mass satisfies M = g n f n at the one-loop level. 
Since this implies g& = 1, in contradiction to experiment, we will allow g s = g n /gA when 
working at the Hartree level. 

In the Hartree approximation, the pion field vanishes for nuclear matter or axially sym- 
metric nuclei, and the energy functional derived from Eq. Q57p reduces to a special case of 
the functional considered earlier. In nuclear matter there are four free parameters: c 2 , c 2 , 
a — a' — rf / g^, and (. We can adapt the analysis of Sec. [JTJ to solve for these parameters 
in terms of the inputs eo, po, K , and Mq. Naively, this would seem to imply that we can 



find good sigma-model descriptions of finite nuclei, in contrast to the conclusion of Ref. |L7 
We first present the analysis and then discuss why this is not the case. 

To specialize the general functional of Sees. || and |T| to the chiral models from Ref. |17 
we set [see Eq. (f?2|)1 

2 

a' = a = ^>0 , (58) 
9t 



3 1 



< , (59) 



3 1 



> . (60) 



Now £ takes the form 



1 "'2 1 ^4 



£ = W p B + E k {<$>- p B ) - —W* - -cw 

+ ^ 2 (l-^ + i ^) + ^ 2 (M 2 -M-). (61) 

At equilibrium, $o still follows trivially from Mq\ 

$ = M - M* , (62) 

and the scalar potential and its derivatives can be evaluated from Eqs. (|26|)-(p8|) using the 
relations ( j59|) and (|60D . We can then simply rewrite the equilibrium conditions from Sec. [1 1 1 
(as before, a subscript "0" denotes the value at equilibrium): 



1. The scalar field equation becomes 

1 



p s0 -aM *iy 2 -i$ 

C s 



3$o , <&§ 1 n M*(M 2 -M* 2 ) 



2M 2M 2 



2c 2 M 2 



(63) 



where p s is still given by Eq. (1321). 
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2. The vector field constraint equation is 



1 £ 



Po - -^Wo + aW [M 2 - Mq 2 } - ^W* = . (64) 



c 2 



3. The energy density can be written explicitly as 

So = (e + M)p 

= \w«p« + + + Po) , (65) 

with ^($o; Po) evaluated from Eq. (|23|). The second equality is obtained from Eq. (|61~D 
using the equation of constraint (|64]). Observe that there is no explicit dependence on 
a in Eq. (|65|) . 

4. The Hugenholtz-van Hove theorem goes through as before, and we again find 

Wo = e + M - \Z(y 2 + I„' 2 • (66) 



5. The compression modulus at equilibrium is 



K = 9p 



71 



_2k F oEp 
with p' s0 still evaluated from Eq. 



v 



M y^ -v / ) ; 



UH-p' s o + l 2 oVo-aW \ 
(H), but now 
Z = 2aM*W , 



(67) 



(68) 



v 



Wo 



Po + K^o 



3 ' 



(69) 



and 



^0 - 2 



$ 3$ 2 - 
M 2M 2 . 



(70) 



One can solve Eqs. d62|)-([70D for c 2 , a, (, and c 2 with some straightforward algebra. After 
eliminating a, (, and c 2 , the resulting equation involving c 2 is actually linear and yields 

2 m 2 + (e; ) 2 ~ e f * ~ ^ U 2 6 + 6e4 1 Wq ^) 
w - 66 j + Wb6 - + ^ 2 

where 
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51 9p 2fc F0 £* ' 
g 

6 = T7f [(e + M)p - £fco] - 3p 
e 3 = Pso + Pso/M* , 

with £ko = £k(po, Mq). The remaining parameters are determined from 




(72) 

(73) 
(74) 
(75) 

(76) 



a 



p s0 - -J $ 



3$ r 



*2 



2.1/ 2M 2 



c = - 



(e + M — - W )po - £fco - ^ $0 



M 4M 2 



p + «W (M 2 - M * 2 ) - C^ 3 /6 



(77) 
(78) 
(79) 



We can use these results to explore the entire range of acceptable normal solutions. The 
procedure to find parameter sets is direct: specify eo, po, Mq, and K Q , and the parameters 
follow algebraically. Nevertheless, the resulting parameter sets will not always be acceptable. 
For example, parametrizations with c 2 < or c 2 < are unphysical. In most cases, however, 
the main issue is the influence of the Lee- Wick solution. 

If we recall from Eq. (|32|) that p s oc M* , it is clear that at any density, Eq. (|63"D will always 
be satisfied when M* = or $ = M. This is the Lee-Wick solution, which is an alternative 
to the "normal" solution, the latter being identified as having zero energy at p B = 0. (In 
general, there is also an "anomalous" solution, which arises from a local maximum in the 
energy density as a function of M* .) There is no evidence that Lee- Wick states exist in 
nature. Notwithstanding, while the Lee-Wick solution exists formally at all densities, we 
are not guaranteed that the normal solution specified at equilibrium will persist at higher 
density (it may coalesce with the "anomalous" solution and disappear). 

An example of the difficulties of bifurcation and coalescence that typically arise from 
the highly nonlinear field equations is found in the original work of Kerman and Miller || . 
In this aoj model [Eq. (]6l|) with a = ( = 0], there are only two free parameters, as in the 
original Walecka model. Proceeding as in that case (see Sec. ||D, one discovers that there 
is no normal solution that reproduces the desired eo and po as an equilibrium point. The 
best one can do is to satisfy one of the conditions, for example, to reproduce the desired 
energy/particle eo at a density po; however, the system will not be in equilibrium (the 
pressure is negative), and moreover, the normal solution fails to exist at densities p > po. 
This situation is depicted in Fig. |5|. Note that the Lee- Wick solution has the lowest energy 
at most densities and is the only solution at high density. 
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FIG. 5. Nuclear matter binding-energy curves for the simple linear auj model. The normal solution has 
the standard eo at density poj but this is not a point of equilibrium. 

These results suggest that we might consider a parametrization acceptable only if it 
produces a normal state with the lowest energy, at least up to the equilibrium density. 
Alternatively, if one recalls that our energy density is truncated at some finite power of 
$/M, one might say that the Lee-Wick solution with $ = M is simply irrelevant. In 
particular, one could argue that it is possible to add (chirally invariant) repulsive terms of 
very high order in <£>/M that would have negligible effect on the normal solution, but which 
increase the energy of the Lee- Wick state enough to raise it above the normal state. Rather 
than digress into a somewhat extraneous discussion about the proper role of the Lee- Wick 
solution in an effective field theory, we will instead impose on our parameter sets a less 
restrictive, more physical constraint that is incontrovertible: If a normal solution is found 
that satisfies the desired inputs, it cannot disappear at too low a density to be useful. 

To make the consequences of this constraint more concrete, so that we can determine the 
boundaries of acceptable parametrizations, we will fix eo and po as usual and then scan the 
(K , Mq) plane, finding the unique normal parameter set (c^, c^, a, and () at each point. 
This allows us to determine regions of the plane in which it is impossible to find a normal 
solution at p > i n which a normal solution exists at p but disappears at some higher density, 
and in which a normal solution exists at all densities. 

The results are summarized in Fig. || where some points corresponding to specific mod- 
els0 from Tables | and |I| are indicated. (See Table |T| for parameter values.) The solid and 



13 These model parameters from Ref. j^] have been adjusted slightly so that nuclear matter 
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FIG. 6. Mq/M vs. Kq for generalizations of the conventional linear sigma model. All results have 
cq = — 16.1MeV and po = 0.1484fm~ 3 . The box denotes the desirable region, and the results of some 
specific parametrizations are indicated by capital letters. The curves are discussed in the text. 

dotted curve is the locus of solutions with ( = 0; on the solid branch, the normal nuclear 
matter solutions exist at all densities, while on the dotted branch, the normal solutions 
disappear at p < 1.38po- (There is indeed a discontinuity in the behavior of the solutions.) 
Regions of the plane to the right and below this curve require ( < 0, which is physically 
unacceptable, as discussed earlier. The dashed curve marks the boundary where normal 
solutions with finite ( remain stable at all densities; thus, only in the small region of the 
plane between the dashed and solid curves can such solutions be found, and this region lies 
far from the desired one. 

The dot-dashed curves are labeled by the value of the density at which the normal 
solution with the specified inputs disappears. Therefore, while it is technically possible to 
find parameter sets for inputs lying within the box, the normal solution disappears at much 
too low a density to be useful. (For example, chiral-model parameters that reproduce the 
inputs at point B generate a normal solution that vanishes at p = 1.07p .) 

Although the choice of limiting density at which the normal state should exist is some- 
what subjective, it is reasonable to demand that it persist at least up to p = 2p , since 
such densities have clearly been created in the laboratory As indicated in the figure, en- 
forcing such a constraint makes it impossible to achieve the desired values of Kq and Mq. 
The underlying problem is the <3? 3 coefficient, with its fixed sign and large magnitude [see 
Eq. (|59|)1, which yields too much attraction. From Fig. ^ we see that K < solutions (with 



saturates at cq = —16.1 MeV. 
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TABLE III. Parameters for the models from Ref. [|l7]] considered in the figures, given in the 
form of dimensionless ratios as in Eq. (|83|). All parameter values have been multiplied by 10 3 and 
Kq is in MeV. The coefficients have been adjusted slightly so that eo = — 16.1 MeV. All models 
have C = 0, and a = a' . 



Model 


1 

2c 2 M' 2 


1 

2c 2 r M 2 


K 

6M 


A 

TA 


o 


M */M 


K 


A 


1.389 


1.815 


0.0 


0.0 


0.0 


0.54 


560 


B 


1.478 


2.309 


0.940 


-0.900 


0.0 


0.63 


230 


E,F 


3.499 


7.315 


-3.499 


0.875 


14.627 


0.78 


680 


I 


2.781 


4.654 


-2.781 


0.695 


7.061 


0.68 


870 


J 


2.470 


3.840 


-2.480 


0.618 


5.406 


0.59 


920 



a = ( = 0) correspond to values of K and Mq far from the phenomenologically desirable 
values. Allowing a = a' > and ( > can produce normal solutions with the desired 
equilibrium properties, but the consequent nonlinearities are too extreme, and the normal 
solution disappears at too low a density. 

For completeness, we also consider the effect of the one-baryon-loop zero-point energy. 
This is strictly applicable only with a renormalizable subset of the lagrangians (namely, 
those with a = ( = 0), but it has been applied more widely in actual calculations pT| , [T7| . 
The change to the energy density is the addition of 

A£ c = ^$ 3 - -^$ 4 + A£ vac (M*) , (80) 
Stt z Sit 2 



where A£ vac (M*) is the one-loop energy defined in Ref. [20|, which starts at 0($ 5 ) and is 
numerically unimportant here. Historically, the motivation for adding AS C was to include 
some repulsion oc $ 3 to cancel the strong attraction in the linear sigma model. It is now 
possible to find a saturating solution, which exists at all densities (see Fig. 59 in Ref. @), 
but it is far from the favored region, since Mq/M « 0.9 is much too large. Allowing 
additional parameters does not change this result. 

Of course, the real test of any parameter set comes in its predictions for the properties of 



finite nuclei. We solve the equations for the finite system using conventional methods |19| , |22 
For chiral models with a linear realization of the symmetry, we require that g s = g n w 13.4 
OI Q 9s — 9-k/ 9a with ~ 1.26. Thus, once we establish c s from the nuclear matter analysis, 
m s is determined (recall that c s = g s /m s ). As noted in Ref. |17| , this has dire consequences, 



since m s is always found to be over 600 MeV, while phenomenologically successful models 
typically have m s « 500 to 550 MeV. Figures [7| and |^ show the effects of varying the scalar 
mass from 500 MeV to 835 MeV with the nuclear matter properties held fixed (models E 
and F); we relax the preceding constraint on g s to obtain the desired m s . The vector mass 



14 We allow the second possibility to compensate for the tree-level value of qa = 1 within the 
model. 
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FIG. 7. Charge densities for chiral model E or F from Table III with three different scalar masses 
are compared to the density extracted from experiment. For comparison, model Tl from Ref. fl36| , which 
provides a good fit to properties of finite nuclei, is also shown. 
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FIG. 8. Charge form factors |F c h(g)| for chiral model E or F from Table [II with three different scalar 
masses are compared to model Tl from Ref. |3(| and experiment. The form factors are simply computed as 
the Fourier transforms of the charge densities in Fig. 0. 
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m v is also held fixed at its experimental value. Large scalar masses produce oscillations in 
the charge density (Fig. [?]) or enhancements in the charge form factor (Fig. ||), which are 
not observed experimentally. Note that the experimental error bars in Fig. |8| are too small 
to include in the figure; the chiral curves are many standard deviations from the data at 
the second maximum. Moreover, the most unfavorable result shown (m s = 835 MeV) is the 
only one with an acceptable coupling, namely, g s ~ g^/gA] using the tree-level value g s = g^ 
produces m s = 1055 MeV. 

Other deficiencies are clearly revealed in calculations of finite nuclei, and this serves to 
validate our choice of desirable nuclear matter properties. Consider, for example, the single- 
particle spectrum near the Fermi surface in heavy nuclei. The large scalar masses required by 
the known value of g s produce significant changes in the effective central potential seen by the 
nucleons. This is illustrated in Fig. |9|. As the scalar mass increases, the slope of the potential 
in the nuclear surface becomes steeper; indeed, if one "averages out" the central oscillations, 
the dotted potential is essentially a square well. Thus, as the scalar mass increases, levels 
with large (nominal) values of the orbital angular momentum i become more deeply bound. 
This is illustrated in Fig. |1(] by models E', E and F, which have scalar masses of 500 MeV, 
650 MeV, and 835 MeV, respectively; the 1/19/2 level, which should be unoccupied, is pushed 
down into the highest filled shell. Moreover, the large value of Mq/M « 0.8 in these models 
reduces the spin-orbit splittings considerably, as can be seen by comparing the lh 9 / 2 -lh u / 2 
splitting in models E and F with the spectra in the successful models B and Tl (from 
Ref. ||36|| ). The conclusion is that the large values of m s and Mq/M obtained in the chiral 



models preclude a correct description of the shell closures in heavy nuclei. 

If we are desperate to obtain a reasonable model of finite nuclei, we might abandon the 
connection to pion physics by decoupling g s from g w , so that m s can be set near 500 MeV. 
This will produce reasonable charge densities and effective central potentials, as verified by 
Figs. |7| and Nevertheless, the models are still restricted to the limited accessible region 
in Fig. H and the consequences are severe. The best that can be done is to reduce Mq/M 
as much as possible, so that the spin-orbit splittings are improved, and the single-particle 
level structure is reasonable. This is illustrated in Fig. |10] by model I (which exists at all 
densities) and model J (which disappears at p « 5p /4); the lower values of m s and Mq 
produce an acceptable level ordering and shell closure. 

Unfortunately, the tradeoff in these models is also clear: the compression modulus Kq 
gets very large, and this disrupts the binding-energy systematics. This is illustrated in 
Fig. ITT], where the deviation in the calculated surface energy (5a 2 ) is plotted as a function 



of the compression modulus for different models. Here 5a 2 is determined by fitting the 
difference in the calculated^ and experimental binding energies of 16 0, 40 Ca, and 208 Pb to 
a form 5E = (8a%)B + (5a 2 )B 2 ^ 3 . (We find 8a\ < 1 MeV for all models. Other analyses give 
qualitatively similar results.) Apparently, models with compression moduli in the desirable 
range yield accurate surface energies, while the chiral models I and J clearly yield surface 
energies that are too large, roughly 40% larger than the observed value (a 2 ~ 18 MeV). 
Thus, by mapping out the accessible region in parameter space and by characterizing 
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A center-of-mass correction is applied in determining the calculated binding energies [23|. 
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FIG. 9. Schrodinger-equivalent central potentials in 208 Pb for model Tl from Ref. |36| and a chiral 
model (E,F) with three scalar masses m s = 835 MeV, 650 MeV, and 500 MeV. 
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and param- 
eters for the others are given in Table [II and in the text. (The scalar mass m s = 500 MeV in models I and 
J.) 



34 



10 



> 

CD 



0J 
<<3 



6 - 



™ (B) 



Walecka (A) 

o 



200 400 600 

K (MeV) 



(J) 



o 



(I) 



o — 



800 



1000 



FIG. 11. Deviations in the surface energy 5a,2 (see text) vs. the compression modulus for the models 



Tl, T2, and T3 from Ref. [B6 and selected models from Tabic [II 



the dependence of the nuclear observables on Mq, K , and m s , we conclude that models 
built upon the conventional linear sigma model and which feature a "Mexican-hat" poten- 
tial cannot reproduce basic nuclear ground-state phenomenology at the Hartree level. The 
nuclear matter analysis of the simplest version (with a = a' = ( = 0) shows that the 
constraints on the scalar self-couplings forced by the form of the scalar potential are incom- 
patible with desirable equilibrium properties. Attempts to generate better properties by 
allowing a = a' > and £ > inevitably lead to Lee- Wick solutions with lower energy and 
to the absence of a normal solution at densities not far above equilibrium. The large value 
of g s forced by the chiral constraint is an additional problem, since it requires the scalar 
mass to be too large. 

A different realization of the chiral symmetry seems more compatible with observed 
nuclear properties and with successful relativistic mean-field models. Here one takes the 
chiral scalar mass to be large to eliminate the unphysical scalar nonlinearities and then 
generates the mid-range attractive force between nucleons dynamically through correlated 



two-pion exchange p2|- |6q1 . In principle, this approach can be realized within the models 
discussed above, albeit at the expense of much greater complexity, since one must first 
construct a boson-exchange kernel containing correlated two-pion exchange and then allow 
this kernel to act to all orders (for example, in a ladder approximation) to determine the 
NN interaction and the resulting nuclear matter energy density |[39|| . This calculation is 
further complicated by the necessity of maintaining chiral symmetry at finite density, which 
is difficult to do when one uses the non- derivative ( "pseudoscalar" ) 7rN coupling implied by 
the linear realization of the symmetry P J57|Jo^ ]. 

Alternatively, this approach can be implemented more easily with a nonlinear realization 
of the chiral symmetry. If desired, the mass of the heavy chiral scalar can be kept finite, 
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so that it plays the role of a regulator that maintains the renormalizability of the model 
||, or the mass can be taken to infinity, so that the chiral scalar field decouples, resulting 
in the effective nonlinear model of Weinberg p2[. The strong scalar-isoscalar two-pion 



exchange can be simulated by adding a low-mass, "effective" scalar field coupled directly 
to the nucleon, and scalar self-interactions can be added to include a density dependence 
in the mid-range NN attractive force. The chiral symmetry is left intact, because the light 
scalar is an isoscalar, and the transformation rules on the nucleon field are nonlinear pS . 
Since the pion mean field vanishes, the resulting chiral mean-field theory produces an energy 
functional just like the ones considered in Sec. [TT]. Moreover, although the coupling strength 
g s of the light scalar is comparable to (as verified by explicit calculation of the correlated 
two-pion exchange [^,|6|]), we no longer require g s = g n and are free to adjust g s within 
a reasonable range. In short, the general energy functionals studied earlier (including the 
Walecka model) are consistent with the underlying chiral dynamics of QCD, although its 
realization in nuclear physics is subtle. In the next section, we illustrate this scenario with 
a specific example from Ref. |36 



V. CHIRAL MODELS THAT WORK 

Here we discuss a chiral model that can successfully describe finite nuclei at the Hartree 
level. As we have seen, variants of the linear sigma model, with spontaneous symmetry 
breaking from a "Mexican-hat" potential, impose strong dynamical assumptions on the 
interactions of the light scalar.^ In contrast, by working with a nonlinear representation 
of chiral symmetry, we remove these constraints, because the chiral-singlet scalar (which is 



called a' by Weinberg |32|) can be decoupled from the other fields. We can then introduce 
a new chiral scalar to incorporate the dynamics in the scalar-isoscalar sector, as discussed 
at the end of the last section. 



We consider the lagrangian from Ref. |36 



C(x) = N(i'fVp + <Mr"75<V - M + g s <f> + ■ ■ )n - -{d»V v - d v V^ 
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(81) 



where qa ~ 1-26 is the axial coupling constant, T>^ = + iv^ + ig^V^ is a chirally covariant 
derivative, and U, v^, and a M depend on the pion field (see Ref. |36|] for details). Here we 
write the nucleon field as iV to distinguish it from the field ip used previously, because N 
transforms nonlinearly under the chiral symmetry. A novel feature of this model is that 
the scale dimension d of the new scalar field S(x) is allowed to differ from unity. The 
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A linear sigma model variation, with spontaneous symmetry breaking from a logarithmic po- 
tential, can provide a good description of finite nuclei, as shown in Ref. [69]. 



36 



scalar fluctuation field is related to S by S(x) = So — 0(x). The form of the lagrangian 
is motivated by requiring that the model satisfy the low-energy theorems of broken scale 
invariance in QCD at the tree level in the effective scalar field. See Ref. [^6j for more 
discussion of the motivation and consequences of Eq. (j8ip. 

This lagrangian illustrates one way to specify vacuum contributions, which at one-baryon- 
loop order modify all powers of 0. The polynomial terms in must be combined with 
corresponding counterterms; in this way, the vacuum contributions are absorbed into the 
renormalization of the scalar polynomial. In principle, there are an infinite number of 
unknown parameters. However, if one insists that the low-energy theorems are satisfied at 
tree level in the meson fields, the end result for the scalar potential must be of the form 
in Eq. (|8~l"D , where the couplings are renormalized 1^6] .0 This potential can be expanded 
as a polynomial in and all coefficients are determined by the three parameters So, d, 
and H q = ?7igCf 2 iS'g/4. One never has to explicitly calculate any counterterms; the nucleon- 
loop effects are automatically included in the scalar potential contained in Eq. fl8"T|) . The 
assumption of naturalness provides an alternative (but not inconsistent) justification for 
limiting the impact of the quantum vacuum at ordinary density to the renormalization of a 
few parameters, as discussed below. 

Thus we have only a few unknown renormalized constants (parameters); in Ref. EH 
these were determined by fitting directly to finite nuclei, using a \ 2 minimization algorithm. 
Although we will not reproduce the results here, the charge densities, single-particle spectra, 
and binding-energy systematics are in good agreement with experiment. Moreover, the 
predicted nuclear matter equilibrium properties (eo, Po, Mq, and K ) fall within the desired 
ranges. 

To connect this model with the present analysis, we present the energy density £ for 
uniform nuclear matter in the Hartree approximation: 

mVo; Ps ) = \mlSld*{(l - ^H-Wl \} + 9,PnVo - icGfcVo)* 



+ V^) m 2 v V 2 + S^; p B ) , (82) 



The parameter sets Tl, T2, and T3 determined in Ref. |36| exhibit naturalness (see Sect. VI), 
as defined earlier; thus, if we expand the scalar potential in powers of 0, terms beyond O(0 4 ) 
are unimportant (although not negligible) at nuclear densities and below, where $/M ^ 
0.4. With this truncation, the energy density becomes a special case of the general form 
considered earlier, and expressions for k, A, and a are given in Table |I[ 

This example illustrates that the general functional in Eq. ([I]) is consistent with the 
spontaneous breaking of chiral symmetry in QCD in a nonlinear realization. The key point 
is that chiral symmetry imposes no constraints on the scalar potential or the scalar-nucleon 
interaction, if a chiral- singlet scalar field generates the midrange attraction. (Of course, 



17 Another way to include vacuum effects is to insist on the renormalizability of the lagrangian. 
Then the first few powers of are fixed by a set of renormalization conditions, and the higher 
powers contain finite, calculable coefficients. 
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all interaction terms must be isoscalars, and the nucleons must obey a nonlinear chiral 
transformation [04 .) The lagrangian in Eq. (|8l|) provides an explicit example of a model 
that is consistent with chiral symmetry (and broken scale invariance) and that produces 
accurate results for finite nuclei in the Hartree approximation; moreover, the coefficients 
are natural, thus allowing for a truncation of the energy functional. In contrast to linear 



chiral models with "Mexican-hat" potentials, the logarithmic potential in Eq. (|82|) allows 
for relatively weak scalar nonlinearities and the dominance of the cubic and quartic terms 
with reasonable values of the scaling dimension d; the coefficients of the scalar terms are 
natural as a result. Although the logarithmic form was motivated by considering the broken 
scale invariance of QCD, we could equally well impose as a phenomenological principle that 
the coefficients in the energy functional be natural, since the subsequent truncation leads 
to similar results. In short, the energy functionals considered in Sec. [TI| are consistent with 
spontaneously broken chiral symmetry, the broken scale invariance of QCD, and the basic 
observed properties of finite nuclei. 



VI. DISCUSSION 

Descriptions of nuclear matter and finite nuclei, which are ultimately governed by the 
physics of low-energy QCD, are efficiently formulated using low-energy degrees of freedom, 
namely, the hadrons. The application of effective field theory methods to nuclear physics 
problems, however, is in its adolescence, with new developments occurring rapidly. In this 
work, we try to build upon the success of relativistic mean-field phenomenology by using 
ideas from density functional theory and effective field theory. One of the key issues is 
naturalness. 

As noted earlier, we have excluded many terms from our energy functional: higher-order 
polynomials in the vector fields and mixed scalar-vector terms, derivative terms, and so on. 
In retrospect, were we justified in neglecting them? An important assumption in applying 
effective field theories, such as chiral perturbation theory, is that the coefficients of terms in 
the lagrangian are "natural", i.e., of order unity, when written in appropriate dimensionless 
units. This assumption permits the organization of terms through a power-counting scheme, 
because one can systematically truncate the expansion when working to a desired level of 
accuracy f3"5fl . We propose an analogous concept of naturalness for the energy functional, 



which will justify the neglect of higher derivatives and powers of the fields when applying 
Eq. (|I]) to nuclei. This is the basic idea that underlies our effective field theory approach to 
the nuclear many-body problem. 

In performing a mean-field analysis P^f4D[] , one can easily identify dimensionless ratios 



that set the scale of individual contributions to the energy. For example, one can rewrite 
the scaled energy density of nuclear matter, £/M A , in terms of the dimensionless ratios 
g v Vo/M = W/M and g s (p/M = $/M, which then become our finite-density expansion 
parameters. If one expresses the nuclear matter energy density in terms of the scaled field 
variables written above [see Eqs. (|20|) — (p2|)] , one finds that the ratios 



2c 2 M 2 ' 2c 2 M 2 ' 6M ' 24 ' 24 



k A C / 

a, ■■■ (83) 
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TABLE IV. Dimensionless ratios of Eq. (p3|) evaluated in the models of Ref. [36]. The absolute 
values are shown, and all have been multiplied by 10 3 . 



Model 


1 

2ciM' 2 


2£M' 2 


K 

6M 


A 

24 


24 


o 


Tl 


1.48 


2.25 


0.02 


0.03 


1.68 


1.16 


T2 


1.65 


2.52 


0.35 


0.12 


1.43 


1.77 


T3 


1.34 


1.95 


0.32 


0.13 


1.44 


0.31 



should all be of roughly equal size (or at least of the same order of magnitude) for our 
expansion to be "natural" . 

In Table |V|, we illustrate these ideas using parameter sets Tl, T2, and T3 from Ref. |36 



These sets accurately reproduce nuclear properties, and we observe that the parameters are 
consistent with our definition of naturalness. The nonlinear scalar parameters tend to be 
somewhat small, but it is clear that one can shift strength between terms with equal powers of 
the meson fields, which follows because both $ and W are roughly proportional to p B . Thus, 
to make concrete conclusions about the importance of nonlinearities, one should include all 
terms consistent with the level of truncation. (Note, however, that the quartic parameter a ' 
was set to zero in Ref. [36]; when it is included in a fit, the optimal parameter set becomes 
more natural.) Observe also the important result that the highest-order terms retained (£, 
A) do not dominate the energy. 

The usefulness of naturalness extends to gradient terms, which do not contribute in 
uniform nuclear matter. Experience with mean-field models applied to finite nuclei shows 
that the derivatives of the fields are small; one finds that typically, | V$|/M 2 and |VW|/M 2 
are (0.2) 2 or less. If we assume mean-field dominance, such that fluctuations around the 
mean fields are small, and the naturalness of the coefficients in the derivative expansion, we 
can truncate the derivative terms at some tractable order. In the models studied here we 
consider only quadratic terms in gradients of the fields. 

Naturalness implies that we should not expect qualitative changes in our description from 
additional higher-order terms. If the physics is consistently dictated by the last terms added, 
the system (or framework) is not natural. Further support for the naturalness assumption 
in our framework comes from extending the model in Ref. |36] to include ^V^V 1 and 
{V^V^Y terms and then repeating an optimization directly to selected properties of finite 
nuclei. The new fits are very close to the fits obtained without these terms. Furthermore, 
contributions to the energy from the new terms are less than 10% of those from the old terms 



at nuclear matter density, and the old coefficients change only slightly in the new fit [40 
Thus contributions from the higher-order terms can be absorbed into slight adjustments of 
the coefficients in Eq. (0). Moreover, since the terms quadratic in gradients of the field are 
small, one expects that adding additional field gradients will produce only minor effects. 

Naturalness also implies that we can exploit the freedom to make nonlinear field redef- 
initions. Although we have no general proof that finite-density observables are unchanged 
under these redefinitions, it is certainly true at the mean-field level {i.e., with classical me- 
son fields) . This freedom allows us to show that our couplings to the nucleon field are more 
general than is apparent at first. For example, if instead of g s ip(t)ip we have g s ip f (<fi)ip , with 
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/(</>) = + • • ., we can redefine the scalar field to recover our standard form. Similarly, we 
can reproduce contact interactions among the nucleons by appropriately choosing nonlinear 
interactions among the meson fields; eliminating the meson fields using the field equations 
then provides the desired products of nucleon densities and currents. Thus there are nu- 
merous ways to write the energy functional that produce equivalent results, provided one 
keeps all orders in the interaction terms. The basic question is which representation leads to 
the most efficient truncation scheme, namely, one with meaningful dimensionless expansion 
parameters, that can be truncated at low order and still provide an accurate description 
of empirical data, using "natural" coefficients. Whereas we cannot say that other ways 



of writing the energy functional will not be useful jTOfl , we have at least found one way 
of parametrizing the interactions (simple Yukawa meson-nucleon couplings with additional 
nonlinear couplings between meson fields) that satisfies the desired constraints. These issues 
will be discussed further in Ref. [§t| . 

One might instead decide to parametrize the nuclear matter energy density by simply 
expanding in powers of the Fermi momentum k F and by adjusting the coefficients in this 
expansion to the desired nuclear matter properties. In practice, this is neither efficient nor 
illuminating. One finds that to reproduce accurately the density dependence found in suc- 
cessful models, one must include many terms in the expansion [roughly through 0(k\ A ) in 
S/p B ]. Moreover, even with this many terms, extrapolations to higher density are problem- 
atic. If all of the coefficients were independent, they would not be sufficiently constrained by 
the few conditions on nuclear matter near equilibrium. In contrast, in models with meson 
fields, the parameters are correlated, and these constraints are effective in restricting the 
relevant region of parameter space. 

If one accepts the validity of our assumption of naturalness, we can compare our trun- 
cated energy functional to one obtained directly from the Hartree (or one-baryon-loop) 
approximation to a given lagrangian. Although both approaches lead to Dirac equations for 
nucleons moving in local meson fields, the parameters in the energy functional and in the 
lagrangian may actually be related in a complicated fashion. The basic premise of density 
functional theory is that it is possible to parametrize exchange and two-nucleon correlation 
effects through local (Kohn-Sham) potentials (or self-energies), and that is how we may 
interpret our local meson fields. 

This approach is useful, however, only for ground-state properties like the total energy 
and particle densities, and in general, the single-particle wave functions and eigenvalues are 
simply mathematical constructs that have no direct relation to observables. Fortunately, rel- 
ativistic calculations including exchange and two-nucleon correlations show that the Hartree 
contributions dominate the self-energies, which are essentially state independent. Thus we 
expect that the large scalar and vector fields seen at the Hartree level will control the 
dynamics, and that the majority of exchange and correlation effects can be included by 
adjusting the nonlinear parameters (and thereby the density dependence) contained in the 
energy functional; this "Hartree dominance" also makes it reasonable to compare the single- 
particle energies with observed spectra, at least for states near the Fermi surface. These 
considerations again indicate the efficiency of writing the energy functional in terms of self- 
interacting meson fields: at normal nuclear densities these fields are small enough compared 
to the nucleon mass to provide useful expansion parameters, yet large enough that exchange 
and correlation corrections can be included, at least approximately, as minor perturbations. 
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Although this mean-field approach allows for an accurate description of many nuclear 
properties and has been applied to various regions of the Periodic Table, it is really just the 
first step in the development of a systematic treatment of the relativistic nuclear many-body 
problem based on effective field theory. There are many limitations and unsolved problems. 
First, there is the question of a more complete inclusion of higher-order corrections in the 
energy functional. These could be computed by starting with a lagrangian whose parameters 
are fit to NN data and then by calculating at a similar level (for example, DBHF) at finite 
density. Although it is known that one can parametrize the density dependence of the 
DBHF self-energies with meson self-interactions [}42[| , and that one can include rearrangement 
effects consistently in the single-particle Dirac equations [^2], the state dependence of the 
self-energies is neglected in these calculations. By comparing results of more sophisticated 
nuclear matter calculations with their mean-field parametrizations, one may be able to 
optimize the description of state-dependent effects, which will be important if one wants to 
study the nucleon-nucleus optical potential at high energy, for example. 

The description of the quantum vacuum should also be improved. As discussed earlier, 
it is possible to include these effects at the one-loop level by incorporating them into the 
model parameters; if our naturalness assumption is valid, only a few such parameters are 
needed in practice. Higher-order contributions, however, will require counterterms involving 
fermion fields and will introduce vacuum contributions that explicitly involve the valence 
nucleons (the strong Lamb shift, for example). Moreover, the present approach neglects 
contributions from meson loops. 

One must also consider the limitations of the present model, particularly with regard 
to extrapolation to high density. At nuclear densities and below, we have seen that only a 
few terms are needed in the energy functional to describe nuclei, and that the sensitivity to 
higher-order products of fields is minimal. This precludes a determination of the coefficients 
of these higher-order terms from empirical data. As the density increases, however, $/M 
and W/M increase accordingly, and the unknown higher-order terms may become important. 
This prevents a controlled extrapolation to high densities, such as those relevant in the inte- 
riors of neutron stars. It is not clear at present what the limits to reliable extrapolation are. 
These difficulties are shared by all existing treatments of the nuclear many-body problem in 
terms of hadrons; for example, calculations of few-nucleon systems based on nonrelativistic 
potentials |7T| , [72| are sensitive to three-body (and perhaps four-body) forces, but reliable 
extrapolation to high density is likely to require information on additional many-body forces 
as well. 

To close the discussion, we return to the question of whether the nuclear properties 
considered here lead to additional constraints on the nuclear matter saturation curve near 
equilibrium. As noted in Sec. |H|, some recent work indicates that the ratio of the "skewness" 
S to the compression modulus K is constrained by monopole vibrations l^jKj. To study 
this issue, we considered an energy functional of the form in Eq. ([[]) with a = a' = and 
varied the ratio (S/K) by varying ( while holding all other nuclear matter inputs fixed. 
We found little correlation between the value of (S/K) and the nuclear observables studied 
here; for example, in Fig. [12|, we show the proton levels near the Fermi surface in 208 Pb 
from five models for which (S/K)o ranges by more than a factor of five. As long as the 
parameter set is natural, the effects of varying the skewness are small, and similar results 
were found for other observables. These results support Bodmer's conjecture [113 that the 
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FIG. 12. Single-particle proton levels near the Fermi surface in 208 Pb for five models with fixed eo, po, 
M*, and if , but different (S/K) . 



spin-orbit splittings in nuclei are determined essentially by Mq and are independent of the 
specific form of the nonlinear interactions. 



VII. SUMMARY 

In this paper, we demonstrate that it is impossible to unite relativistic mean-field phe- 
nomenology with manifest chiral symmetry by building upon the conventional linear sigma 
model, which contains a "Mexican-hat" potential and a light scalar meson playing a dual 
role as the chiral partner of the pion and the mediator of the intermediate-range NN attrac- 
tion. We illustrate the generic failure of this type of model at the Hartree level as well as 
the characteristics of chiral models that can successfully describe finite nuclei. 

Although our study of chiral models involves specific lagrangians treated at the Hartree 
(or one-baryon-loop) level, the analysis is based on the more general concept of an energy 
functional that describes the finite, many-nucleon system. The energy functional depends 
on valence-nucleon Dirac wave functions and classical scalar and vector meson fields. We 
interpret this energy functional within the framework of density functional theory, which says 
that the exact energy functional for the system can be written entirely in terms of the scalar 
density and baryon four-current density. Rather than work solely in terms of densities, we 
include auxiliary meson fields, which correspond to (local) Kohn-Sham potentials, and which 
provide an efficient parametrization of the density dependence. Thus, with a sufficiently 
general energy functional, we can include many-body effects beyond the simple Hartree 
level, even though we retain only classical meson fields and local interactions. Rather than 
attempt to compute the form of the functional directly from an underlying lagrangian, we 
use observed nuclear properties (and extrapolations to nuclear matter) to determine the 
unknown coefficients. 

In the general case, the energy functional contains an infinite number of interaction terms 
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involving the nucleon and meson fields, and their derivatives. This is actually an advantage, 
since transformations of the field variables allow one to cast the functional in different forms, 
with the goal being to find the representation that allows the most efficient truncation. For 
the truncation to be plausible, one must have a reasonable (small) expansion parameter, 
and one must assume a principle of "naturalness", which says that the coefficients of the 
expansion, when written in appropriate dimensionless ratios, must all be of order unity. 
These are precisely the basic ideas underlying applications of effective field theories, such 
as chiral perturbation theory. We showed that such a truncation of the energy functional 
is possible at and below normal nuclear densities, if the functional is written in terms of 
self-interacting scalar and vector fields with Yukawa couplings to the nucleons; under these 
conditions, the ratios of the meson fields (and their gradients) to the nucleon mass are 
small enough to serve as useful expansion parameters, and sets of coefficients can be found 
that both reproduce observed nuclear properties and satisfy the naturalness assumption. In 
practice, truncation of the energy functional at quartic terms in the meson fields and at 
quadratic terms in field gradients is sufficient for the observables studied here. 

To constrain the energy functional, we use nuclear observables that are meaningful from 
the standpoint of the underlying density functional framework, such as (charge) densities and 
total binding energies. Although in general, the single-particle Dirac eigenvalues obtained 
in the local potentials are not directly related to observables, we note that in the relativistic 
nuclear many-body problem, the single-particle potentials (or self-energies) are known to 
be dominated by the Hartree contributions and are essentially state independent. Thus the 
exchange and correlation effects included in the potentials are small, and it is reasonable 
to also use observed nuclear spectra as input, at least for states near the Fermi surface. In 
short, the observables we use to constrain the parameters are precisely the ones that should 
be calculable in a density functional framework where the Hartree terms dominate the self- 
energies. We also observe the fortunate result that the scalar and vector fields are small 
enough compared to the nucleon mass to serve as useful expansion parameters, yet large 
enough that exchange and correlation effects enter only as small perturbations (at least for 
occupied states). 

Rather than fit directly to the properties of finite nuclei, we extrapolate instead and 
determine our parameters from the properties of nuclear matter near equilibrium. In nu- 
clear matter, the energy functional reduces to a function of the Fermi momentum and the 
(constant) scalar and vector fields; it is then possible to invert the field equations and ex- 
press the unknown parameters directly in terms of the nuclear matter "observables". This 
allows us to scan the parameter space to find sets that lead to an accurate reproduction of 
nuclear properties. More importantly, we can exclude models that are too constrained (by 
symmetry, for example) by showing that it is impossible for them to reproduce desirable 
nuclear matter results. 

Parameter sets that accurately reproduce nuclear observables typically involve small but 
significant nonlinear meson self-interactions. These self-interactions introduce additional 
density dependence beyond that contained in simple Yukawa meson-nucleon couplings. This 
density dependence can be interpreted either in terms of meson masses that are modified 
in the nuclear medium or in terms of variations in the scalar and vector self-energies. We 
can make three important observations about this additional density dependence: First, it 
can be parametrized efficiently using self-interactions of the meson fields. Second, to make 
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meaningful statements about the size of these self-interactions, one must include all allowed 
terms at the given level of truncation, as required by naturalness. This is relevant since most 
relativistic mean-field calculations to date have concentrated on scalar self-interactions only, 
omitting scalar-vector and vector-vector terms. Third, although we choose to determine 
these interactions by fitting to empirical nuclear properties, one could instead fit them to 
the results of more detailed nuclear matter calculations based on an underlying lagrangian, 
such as DBHF. If the DBHF self-energies are determined from fits to NN observables, this 
alternative procedure provides a direct link between relativistic finite-nucleus structure cal- 
culations and the two-nucleon problem. 

To address the chiral models, we considered specific lagrangians that realize the sponta- 
neously broken chiral symmetry of QCD in different ways and studied them at the Hartree 
level. The resulting energy functionals are special cases of the general one studied earlier, 
so we apply our previous analysis and test whether acceptable parametrizations can be ob- 
tained. We find that the accessible region for linear sigma models is far from the region of 
acceptable nuclear phenomenology. The results of this new analysis solidify our previous 
conclusion ]T7[ that chiral hadronic models built upon the conventional linear sigma model 
cannot reproduce observed properties of finite nuclei, at least at the Hartree level. The two 
critical constraints imposed by the conventional lagrangian are the "Mexican-hat" potential 
and the connection between the scalar coupling and the pion coupling (g B = g^jg^). The 
former generates scalar self-interactions with the wrong systematics, and the latter produces 
a scalar mass that is too large, leading to unrealistic nuclear charge densities, level orderings, 
and shell closures. 

On the other hand, one can find chiral models that are phenomenologically successful 
at the Hartree level. As an example, we discuss a model that was introduced in a recent 
paper, which features a nonlinear realization of chiral symmetry. One motivation for such 
a model is that the chiral scalar present in the linear sigma model is irrelevant for the 
nuclear dynamics; the mid-range NN attraction arises instead from the exchange of two 
correlated pions. One can simulate this attraction by introducing an auxiliary scalar field 
with a Yukawa coupling to the nucleon. Chiral symmetry is maintained if the nucleon 
field transforms nonlinearly and the auxiliary field is a chiral scalar. Due to the nonlinear 
realization of the symmetry, the "Mexican-hat" potential is no longer needed, and various 
constraints on the model parameters disappear. At the Hartree level, the energy functional 
has the full flexibility of the general functional considered earlier; thus, parameter sets can 
be found that accurately reproduce nuclear properties in a framework that is consistent with 
the underlying chiral symmetry of QCD. Moreover, we observe that successful models have 
parameters that are consistent with our "naturalness" assumption, which forms the basis 
for our expansion and truncation scheme. 

Although the general mean-field approach studied here leads to excellent results for the 
bulk properties of nuclei, there are many questions to be answered before one has a consis- 
tent, systematic treatment of the relativistic nuclear many-body problem. A more accurate 
inclusion of vacuum, exchange, and correlation effects should be pursued. The extension of 
this framework to deal with high-density matter, where the expansion parameters are no 
longer small, is also an outstanding problem. Some of these ideas will be discussed in a 



future paper [40 
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